Investigation of the reaction pp — > pK n + A in search 
of the pentaquark 

M Nekipelov 1 ' 2 , M Biischer 1 , M Hartmann 1 , I Keshelashvili 13 , 
V Kleber 4 , V Koptev 2 , Y Maeda 5 , S Mikirtychiants 2 , R 
Schleichert 1 , A Sibirtsev 6 , H Stroher 1 and Yu Valdau 1 ' 2 

1 Institut fiir Kernphysik, Forschungszentrum Jiilich, 52425 Jiilich, Germany 

2 Petersburg Nuclear Physics Institute, 188300 Gatchina, Russia 

3 High Energy Physics Institute, Tbilisi State University, 0186 Tbilisi, Georgia 

4 Physikalisches Institut, Universitat Bonn, 53115 Bonn, Germany 

5 Research Center for Nuclear Physics, Osaka University, Ibaraki, Osaka 567-0047, 
Japan 

6 Helmholtz-Institut fiir Strahlen- und Kernphysik, Universitat Bonn, 53115 Bonn, 
Germany 

E-mail: m . nekipelov@f z- juelich . de 

Abstract. The reaction pp — » pK°Tr + A has been studied with the ANKE 
spectrometer at COSY- Jiilich at a beam momentum of 3.65GeV/c in order to search 
for a possible signal of the pentaquark O + (1540), decaying into the pK° system. By 
detecting four charged particles in the final state (ir + , ir~ and two protons), the K° 
and the A have been reconstructed to tag strangeness production. It has been found 
that the n + A missing-mass spectrum displays no significant signal expected from the 
O + (1540) excitation. The total cross section for the reaction pp — > pK°ir + A has 
been deduced, as well as an upper limit for the 9 + production cross section. The 
intermediate A ++ K°A state seems to provide a significant contribution to the reaction. 
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1. Introduction 

The reactions pp^NKnA have been intensively used for hadron spectroscopy at beam 
momenta above 5 GeV/c, since they offer an excellent opportunity to search for baryon 
and hyperon resonances by studying 7rN and it A final states, respectively [1-9]. It 
also allows one to investigate the KN system, which attracted a lot of attention due 
to its putative coupling to the pentaquark baryon + (154O) [10]. Despite a number 
of experimental indications for the existence of such a state, recent negative results, 
mostly obtained at higher energies (see e.g. [11] and references therein), have cast 
serious doubt on the existence of the pentaquark © + . Its current status is reviewed 
in [12]. There is no theoretical investigation that reconciles both the positive and the 
negative observations, although it has been shown in [13] that the production of + 
should be strongly suppressed at higher energies. Hadronic experiments at low energies 
are believed to be crucial to clarify the situation. 

The main features of the available data for pp^pK°iT + A [1-9] are as follows: 
a) Structures in the tt + A and ir + p systems have been clearly identified as £+(1385) 
and A ++ (1232) resonances, respectively, providing the dominant contributions to the 
reaction cross section. In the n + K system the K*(892) resonance has been detected, 
although its contribution is rather small, b) An enhancement in the KA invariant mass 
spectrum at 1.7 GeV, which is expected because of the excitation of the Pn(1710) and 
Pi 3 (1720) baryon resonances [12], has been observed [3,4]. In [1] it has been argued that 
this enhancement results from a kinematical reflection due to S + (1385) production, c) 
The KN and A A" invariant mass spectra have not been investigated yet. All published 
data suffer from rather limited statistics and a moderate energy resolution. 

It is thus important to re-examine this reaction with new measurements with 
improved statistical accuracy and high mass resolution. Such a dedicated experiment has 
been performed with the magnetic spectrometer ANKE [14] at the COoler SYnchrotron 
COSY [15] at the Research Centre Julich. 

2. The experiment 

COSY provided a circulating proton beam of p p = 3.65GeV/c with an intensity 
of approximately 4 x 10 10 protons per spill interacting with a hydrogen cluster jet 
target [17] with a thickness of about 5 x 10 14 cm -2 . The detection systems of ANKE, 
capable of detecting charged particles with angles 6 < 12°, were triggered by a three- 
fold coincidence of two positively charged particles and a n~ (from A decay). The 
reaction has been identified by the off-line selection of four particles, namely a tt + and 
7r~ meson and two protons. Detection and identification of all the charged particles 
was performed by a combination of scintillation counters and multi-wire proportional 
chambers (MWPC). While the first allowed one to exploit the time-of- flight technique, 
the MWPCs supplied full tracking and momentum information, and permitted the 
selection of particles originating from the target via analysis of their vertical angle (for 
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Figure 1. a) Mass contour plot: invariant mass of ir~p versus missing mass of 
ir + ir~pp. An enhancement corresponding to the masses of K° (497.7 MeV/c 2 ) and 
A (1115.7 MeV/c 2 ) is seen. Corresponding projections are presented in figures b and 
c. A lego plot is given in figure d. Shaded areas in the projection plots onto the axes 
ininv(^ p) (c) and m x (ir + ir~ pp) (b) correspond to the independent cuts on masses 
of K° and A, respectively. The dashed lines indicate these cuts. The hatched areas, 
marked as 1 and 2 in the mi nv (ir~ p) plot (c) denote the regions used for the background 
subtraction (see text). 

details see [14,16]). A coincidence measurement of four particles reduced the accidental 
background to a negligible level. The amount of misidentified particles was also small, 
due to the analysis of time-of-flight differences between detected particles. Therefore, 
the remaining background was almost entirely of physical origin. 

The reaction pp —> pK°ir + A has been reconstructed by identifying the A hyperon 
in the ir~p invariant mass and the K° meson in the missing mass of the ir + 'K~pp system. 
The K° signal thus includes both contributions from K§ and K®. The corresponding 
invariant mass versus missing-mass scatter plot is shown in figure Hk. An enhancement 
of K°A coincidences is clearly visible on top of a sizeable background, remaining after 
cuts on masses are made (see shaded areas in figures [£b and c) . This background was 
removed by the side-band subtraction method. For the subtraction, events in the regions 
besides the A peak are selected, assuming that the background behaviour is smooth. The 
latter assumption was checked by choosing two different background regions (hatched 
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Figure 2. Side band subtraction technique, a) The difference of two missing 
mass spectra is plotted, in which different background regions were selected for the 
subtraction, b) K° signal after background subtraction has been applied. The dashed 
line shows the result of Monte Carlo simulations assuming a pK a ir + A phase space 
behaviour. 

areas marked as (1) and (2) in figure Wp), and building the bin-by-bin difference between 
them. The side-band subtraction method can be applied, if such a difference spectrum 
does not exhibit any distinctive structures. For our case the result of this investigation 
is presented in figure Eh for the missing mass distribution m(7r + A), indicating that no 
signal shape distortion should occur because of background subtraction. The resulting 
K° missing mass spectrum after background subtraction is shown in figure [2b. The 
dashed line shows the result of Monte Carlo simulations assuming a pK rr + A phase 
space behaviour, which reproduces the experimental spectrum. 

Differential distributions of the individual particles, such as momenta and scattering 
angles, are also well described by Monte Carlo simulations with the total reaction cross 
section as a free parameter and assuming four-body phase space distributions in the 
final state. The missing mass distribution m(n + A), shown in figure [31 is the one where 
the signal from the pentaquark would be expected to appear. However, no obvious 
peak at 1.54GeV/c 2 is visible. The total number of collected pK°n + A events is 1041, 
which corresponds to a rate of 82 events/ day for a total integrated luminosity of 36 pb _1 . 
Together with the experimental data, results from the simulations are plotted as the 
solid line in figure [3) As one can see, only a moderate agreement (xidf = 1-41) between 
experimental data and simulations based on four-body phase space is achieved. While 
for the high mass part of the spectrum the agreement is reasonable, an obvious excess 
of experimental events is observed around 1.47GeV/c 2 . Such an enhancement may be 
connected with the excitation of some intermediate resonances. 

Indeed, at higher energies the reaction pp —>■ pK°n + A proceeds mainly through 
the production of intermediate resonances A ++ (1232) and S + (1385) [2,4]. It has also 
been argued that A r *(1650) or AT*(1710), as well as an excitation of if* (892), might be 
important [4]. 
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Figure 3. Missing mass spectrum of ir + A system. The solid line shows the result of 
Monte Carlo simulations assuming pK°ir + A phase space behaviour. 



In order to check these hypotheses a simultaneous fit of the two mass distributions, 
namely m x (7r + A) on one hand and m llw {TT + p) (for A ++ (1232)), or m inv (7r + A) (for 
E + (1385)), or m x (7r + p) (for iV*(1710)) on the other, has been performed. Since the beam 
momentum is close to the K*(892) production threshold the relevant mass region is not 
populated, and therefore this resonance does not contribute significantly. By introducing 
the TV* (1710) resonance into the simulations the agreement with experimental data 
significantly decreases, thus unfavouring this channel. No significant difference is found 
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Figure 4. a) Invariant mass spectrum of 7r + A system. The dashed lines restrict the 
area where a signal of £(1385) is expected. 



in the quality of the fit including both the non-resonant channel and the one with an 
intermediate S + (1385), although a small amount of events can be attributed to the 
production of S + (1385) (see figure Hj). The fit results in a S + (1385) peak area, which 
is compatible with being a statistical fluctuation of the background. Therefore, only an 
upper limit for the possible S + (1385) production can be deduced. Our estimates are 
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based on the method of calculating the upper limit for the peak area as given in [18]. 
This method assumes that the maximum peak area depends on the background and the 
total counts in the peak region. The upper limit is then obtained for a given significance 
level, taking into account that the background is not exactly known and has a certain 
error. As a result we extract N% < 69 events at the 95% confidence level, from which 
an upper limit for the cross section is deduced, given further below. 

The contribution from the formation of an intermediate A ++ (1232) becomes 
obvious in the mass distribution of the 7r + p system shown in figure [5k.. Here the dotted 
region originates from the pK°n + A phase space, while the solid line depicts the sum of 
the latter and contributions from A ++ (1232) and G + . The overall agreement between 
experimental data and simulations improves (Xndf = 1-24), as shown in figure [5b. The 
inset in figure [5b presents the individual contributions from non-resonant production 
and from intermediate A ++ (1232) excitation. Since the shapes of the distributions above 
1.47GeV/c 2 are very similar, any conclusion about the possible 6 + production does not 
depend on the assumption made for "non-exotic" mechanisms. 
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Figure 5. a) Invariant mass of n + p system. The dotted area is obtained from the 
simulations assuming phase space without intermediate resonances, and the solid line 
is the sum of all the contributions, b) Missing mass spectrum of 7r + A system. The solid 
line denotes the sum of all the contributions. The inset shows individual contributions 
from non-resonant production (dotted) and from intermediate A ++ (1232) (hatched). 
The black region corresponds to the maximum permissible + (154O) signal. 



Finally, the black area in (figure [5b) corresponds to the possible signal expected 
from the + production. It has been found that the maximum permissible number of 
pentaquark events weakly depends both on the width and the position of the resonance. 
The width was varied from 1 up to 15MeV/c 2 , with the latter value comparable to the 
missing mass resolution (FWHM) of the apparatus, while the + mass ranged between 
1.52 GeV/c 2 and 1.54 GeV/c 2 . The maximum number of pentaquark events found under 
these various assumptions is given by the best fit to the experimental data and equals 
to 28 ± 20. It is obvious, that even with the best fit the signal is compatible with zero 
within errors, and therefore cannot be used for the cross section estimates. Thus, as in 
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the E + (1385) case, only an upper limit for the cross section will be given below. For the 
evaluation of the upper limit of the + yield at the 95% confidence level the approach 
from [18] has again been used. This resulted in 44 + events, used for calculation of an 
upper limit of the cross section (see below). 

For normalisation of the data, the reaction pp — > pK + A has been used. This 
reaction was measured simultaneously during the experiment. Again all four particles in 
the final state have been detected, including the proton and n~ from A decay. Contrary 
to the pp — > pK°7r + A reaction there are no background channels in the case of the 
production of pK + A (see figure [6]). Such a normalisation procedure, where a reaction 
with known cross section [19, 20] and nearly the same final state is measured, allows 
one to avoid most of the systematic uncertainties connected with particle identification 
and detector efficiencies, although adding a contribution due to the uncertainty in the 
acceptances of both measured channels. 
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Figure 6. Missing mass spectrum of K + p system in the reaction pp — > pK + A. A 
background-free final-state identification has been achieved. 




3. Results and discussion 

Since the acceptances of both four-body phase-space production and production with 
an intermediate A ++ (1232) are very similar, a total cross section for the pK 7r + A final 
state can be calculated largely independent of the decomposition into separate channels. 
After the normalisation and efficiency corrections, the following total cross section for 
the pK°n + A final state has been deduced: 

a tot = 1.41 ±0.05 ±0.33 /ib, 

where the first error is statistical, while the second is systematic. The systematic 
uncertainty is mostly coming from the error of <j p k+a [19]. 

The total cross sections for the non-resonant channel and the channel with the 
A ++ (1232) excitation have also been evaluated: 

(T p n K0~lT nant = °- 92 ± 0.16 ± 0.21 fjb, 
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a A++ ^o A = 0.49 ± 0.14 ± 0.11 fib. 

Unlike with the total cross section, where the calculation relies on the number of 
experimental events, the results for the individual channels include an additional error 
from the decomposition. The measured cross section for the pp — > A ++ K°A reaction 
is significantly lower than a model prediction, a ~ 6 yub [20] . However, this model 
overestimates as well the data available at high energies. 

Table 1. Contributions of different channels to the cross section for pp — > pK°ir + A 
in per cent of the total cross section cr t ot (see also figure [7]) . NR denotes non-resonant 
production. 
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Figure 7. Total cross sections for the reaction pp — > pi4T°7r + A as a function of the 
beam momentum. The ANKE data point is denoted by a circle, while the other cross 
sections (see [1-9]) are given by squares. The dashed line corresponds to four-body 
phase space normalised to the ANKE data point. 

The energy dependence of the total cross section of the reaction pp — > pK°n + A 
is shown in figure [71 The data point from the current work is denoted by the full 
circle, while cross sections from other experiments are shown as squares. A nearly 
constant cross section above p p ~ 6GeV/c can be understood if the reaction proceeds 
via intermediate resonances. Otherwise, a steep rise of the total cross section due to 
four-body phase space should be expected, as is shown by the dashed line in figure [7J 
Such a resonance dominance has been observed at higher beam energies [2,4], where 
A ++ (1232) and S + (1385) provide the dominant contribution to the yield. Nevertheless, 
below p p ~ 6 GeV/c, where the energy available for the production of the resonances is 
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limited, non-resonant production is still competitive, and its relative contribution should 
decrease with increasing energy. The opposite behaviour is expected for resonance 
channels. Their contribution should increase with the beam momentum as illustrated 
in table [TJ 

Taking into account the estimated maximum number of 44 events for the G + from 
the discussion above, we have evaluated an upper limit for the cross section of + - 
production in pp — ► + 7r + A at p p = 3.65 GeV/c: 

ve+n+A < 0.058 fib (95% CL) . 

It is not possible to directly compare this result with the cross section for pp — > pi^°E + 
measured at p p = 2.95GeV/c by the COSY-TOF collaboration [21], indicating a + - 
signal. The exit channels are distinct and even the number of particles in the final state 
differs. As a crude estimate, one may assume that production of an additional pion 
leads to the reduction of the cross section by an order of magnitude. However, a similar 
factor can be regained due to the production of a A instead of a E. But a difference 
in beam momentum of 700MeV/c makes such a comparison very rather simplistic and 
questionable. 

The upper limit for the E + (1385) production cross section is also provided by our 
data: 

^s+(i385)^o P < 0.15/ib (95% CL). 

As has already been mentioned above, the estimates of these cross sections do not 
depend on how the corresponding mass distributions are decomposed, since the shapes 
of the separate contributions are very similar in regions where the peaks from + and 
E + (1385) are expected. 

In summary, we have presented measurements of the reaction pp — > pK°n + A for 
p p = 3.65GeV/c. While at higher beam momenta the formation of this final state 
is dominated by the production of intermediate resonances, non-resonant production 
seems to prevail at COSY energies, although the intermediate (A ++ i^°A)-state provides 
an important contribution to the yield. Unlike at higher energies, the contribution of 
E + (1385) is small. We have found no obvious signal indicative of the + (154O). Its 
cross section in pp — > + tt + A reaction is limited to be less than 58 nb. 

Acknowledgments 

We wish to acknowledge the assistance we received from the COSY and ANKE crews 
when performing these measurements. This work has partially been supported by 
BMBF, DFG, Russian Academy of Sciences, and COSY FFE. 

References 

[1] Alexander G et al 1967 Phys. Rev. 154 1284 

[2] Bierman E, Colleraine A P and Nauenberg U 1966 Phys. Rev. 147 922 



Investigation of the reaction pp — > pK°7Y + A in search of the pentaquark 10 

[3] Klein S et al 1970 Phys. Rev. D 1 3019 

[4] Firebaugh M et al 1968 Phys. Rev. 172 1354 

[5] Alexander G et al 1968 Nuovo Cimento 53 A 455 

[6] Klein S University of California Lawrence Radiation Laboratory Report UCRL 18072 (unpublished) 

[7] Dunwoodie W M 1968 University of California Los Angeles UCLA-1033 (unpublished) 

[8] Holmgren S et al 1967 Nuovo Cimento 51 A 305 

[9] Caso C et al 1968 Nuovo Cimento 55 A 66 
[10] Diakonov D, Petrov V and Polyakov M 1997 Z. Phys. A 359 305 
[11] Battaglieri M et al 2006 Phys. Rev. Lett. 96 042001 

[12] Yao W-M et al (Particle Data Group) 2006 J. Phys. G: Nucl. Part. Phys. 33 1 

[13] Titov Al et al 2004 Phys. Rev. C 70 042202 

[14] Barsov S et al 2001 Nucl. Instrum. Methods A 462 364 

[15] Maier R 1997 Nucl. Instrum. Methods A 390 1 

[16] Biischer M et al 2002 Nucl. Instrum. Methods A 481 378 

[17] Khoukaz A et al 1999 Eur. Phys. J. D 5 275 

[18] Hclene O 1983 Nucl. Instrum. Methods A 212 319 

[19] Louttit R I et al 1961 Phys. Rev. 123 1465 

[20] Tsushima K et al 1999 Phys. Rev. C 59 369 

[21] Abdel-Bary M et al 2004 Phys. Lett. B 595 127 



